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ABSTRACT 

We report on the observation of the bright, long gamma-ray burst, 
GRB 090926A, by the Gamma-ray Burst Monitor (GBM) and Large Area Tele- 
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scope (LAT) instruments on board the Fermi Gamma-ray Space Telescope. 
GRB 090926A shares several features with other bright LAT bursts. In par- 
ticular, it clearly shows a short spike in the light curve that is present in all 
detectors that see the burst, and this in turn suggests that there is a common 
region of emission across the entire Fermi energy range. In addition, while a 
separate high-energy power-law component has already been observed in other 
GRBs, here we report for the first time the detection with good significance 
of a high-energy spectral break (or cutoff) in this power-law component around 
1.4 GeV in the time- integrated spectrum. If the spectral break is caused by opac- 
ity to electron-positron pair production within the source, then this observation 
allows us to compute the bulk Lorentz factor for the outflow, rather than a lower 
limit. 

Subject headings: gamma rays: bursts 



Introduction 



Gamma-Ray Bursts (GRBs) are the most energetic transients in the universe. The first 
brief and intense flash, the so-called prompt emission, has been observed in the X-ray and 
gamma-ray bands, while subsequent long-lived afterglow emission has so far been observed 
mainly at energies in the X-ray band and below. The prompt emission is thought to be 
produced in an ultra-relativistic outflow, but its detailed emission mechanism has been a long- 
standing problem. It has been widely believed that the afterglow is the synchrotron emission 
from the forward shock that propagates in the external medium, but Swift observations 
have pointed out some difficulties in this model (for recent reviews, IZhangl l2007t iMeszaros 
20061 ). The study of the gamma-ray emission in the GeV energy range is expected to give 
us important information on these issu es and even on the nature of the progenitors and 
the u ltra-relativistic outflows of GRBs ( IBand et al.ll2009t Falcone et al.ll2008t Fan &: Piran 
2008h . 



The Fermi Gamma-ray Space Telescop e hosts two instrum ents, the Large Area Tele- 
scope (LAT, 20 MeV to more th an 300 GeV; lAtwood et al.ll2009f) and the Gamma-ray Burst 



Monitor (GBM, 8 keV-40 MeV; Meeean et al. 



20091 ). which together are capable of measur- 



ing the spectral parameters of GRBs across seven decades in energy. Since the start of science 
operations in early August 2008, the Fermi LAT has significantly detected 16 GRBs. These 
events, including the very bright long-duration and short-duration bursts GRB 080825C, 
GRB 080916C, GRB 081024B, GRB 090510, and GRB 090 902B, have reveale d ma ny im- 



portant, seemingly common, features of GRB GeV emission (lAbdo et al.ll2009al Jb 



2010 
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Ackermann et al.ll2010al Jbl): (1) the GeV emission onsets of many LAT GRBs are delayed 
with respect to the MeV emission onsets; (2) some LAT GR Bs have extra hard components 
apart from the canonical Band function (IBand et al.lll993l ). which typically peaks in vF v 
between around 100 keV-1 MeV; (3) the GeV emission lasts longer than the prompt MeV 
emission, showing power-law temporal decays at late times. 

In this paper, we report on the analysis of the bright, long GRB 090926 A detected by 
Fermi LAT/GBM. The light curve of this burst above 100 MeV shows a sharp spike with a 
width of 0.15 s, fast variability that we use to constrain the origin of the high-energy photons 
within the spike. Furthermore, from the detection of a break in the > 100 MeV gamma-ray 
spectrum, we derive constraints on the bulk Lorentz factor and the distance of the emitting 
region from the central source. Section [5] summarizes the detections of GRB 090926A by the 
GBM and the LAT, and the follow-up observations. Section [3] presents the light curves of the 
prompt emission as seen by both instruments and describes a sharp pulse seen in all detectors. 
In section HJ we detail the spectral analysis of the burst through time-resolved spectroscopy, 
the measurement of a break in the extra-component, and the extended emission found in 
the LAT data out to 4.8 ks after the trigger. These last two points are at the center of the 
physical interpretation of the observations that is developed in section 15.11 Throughout this 
paper, we adopt a Hubble constant of Hq = 72 km s _1 Mpc -1 and cosmological parameters 
of tt A = 0.73 and Q M = 0.27. 



2. Observations 

At 04:20:26.99 (UT) on Sept. 26, 2009 (hereafter T = 275631628.98 s mission elapsed 
time), the Fermi Gamma-ray Burst Monitor (GBM) triggered on and localized the long 



GRB 090926A at (RA, Dec) = (354.5°, -64.2°), in J2000 coordinates flBissaldill2009l : lUehara et al 



20091 ). This position was ~52° with respect to the LAT boresight at the time of the trigger 
and well within the field of view. An Autonomous Repoint Request (ARR) was gener- 
ated, but the spacecraft initially remained in survey mode as the Earth avoidance angle 
condition was not satisfied by the burst pointing direction. The on-board GBM position of 
GRB 090926A was occulted by the Earth at roughly T + 500 s until it rose above the horizon 
at approximately To + 3000 s. At that time, the spacecraft slewed to GRB 090926A and kept 
it close to the center of the LAT field of view until T + 18000 s, though the source location 
was occulted by the Earth several times over that time period. 

Emission from GRB 090926A was evident in the Fermi LAT raw trigger event rates, and 
the number of LAT events (~ 200 photon candidates above 100 MeV) is comparable to that of 
the other bright LAT bursts, GRB 080916C, GRB 090510 and GRB 090902B. The increase 
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in the photon count rate during the prompt phase is spatially and temporally correlated 
with the GBM emission with high significance, and extended emission is observed until 
T + 4800s. The best LAT on-ground localization is (RA, Dec) = (353.56°, -66.34°), with 
a 90% containment radius of 0.07° (statistical; 68% containment radius: 0.04°, systematic 
error is less than 0.1°) and is consistent with the XRT localization. 

Indeed, based upon the GCN report issued for the LAT detection, a Swift TOO ob- 
servation was performed, and an afterglow for GRB 090926A was detected with XRT and 
UVOT at T + 47 ks and locali zed at (RA, Dec) = (353.40070°, -66.32390°) with an uncer- 
tainty of 1.5 "(90% confidence) (jVetere et al.ll2009f ). VLT observations determined a redshif t 
for GRB 090926A of z = 2.1062, using the X-shooter spectrograph flMalesani etaD l2009h . 
Suzaku/WAM and Skynet /PRO MPT also detecte d the soft gamma-ray prompt and optical 
afterglow emission, respectively (INoda et al.ll2009l ). 



3. Light Curves 



In figure (TJ we show the GBM and LAT light curves in several energy bands. The 
highest energy photon is a 19.6 GeV event, observed at To + 25 s within 0.03° from the LAT 
position of GRB 090926A, well within the 68% containment of the point spread function at 
that energy. The light curves show that the onset of the LAT e mission is delayed b y 3.3 
s with respect to the GBM em ission, similar to other LAT GRBs JAbdo et al.ll2009afl clldl. 



2010; Ackermann et al. 



2010a 



bj). Deta iled analysis of the GBM data results in a formal 



T90 duratiorfl (jKouveliotou et al.lll993l ) of 13.1 ± 0.2 s, with a start time of T + 2.2 s and 



a stop time of T + 15.3 s. The emission measured in the LAT above 100 MeV has a similar 
duration; however, owing to the efficient background rejection applied to the LAT data, the 
signal is clearly visible in the light curve well after this time range. 

The time intervals chosen for spectroscopy are indicated by the vertical lines in figure HJ 
with boundaries at T o +(0, 3.3, 9.8, 10.5, 21.6) s. The end of the last time interval at T +21.6s 
was chosen somewhat arbitrarily as the end of the prompt phase, but we carefully verified 
that our results are not affected by a slightly different choice. Figure |2] shows a zoom of 
some of the light curves between T + 2.2 s and T + 15.3 s, with a binning of 0.05 s, and 
highlights the presence of the sharp peak seen in each of the Nal, BGO, and LAT light 
curves at T + 10 s. As seen on figure [2J the peak is clearly in coincidence in all of the light 
curves, indicating a strong correlation of the emission from a few keV to energies > 100 



1 The T90 duration is the time over which the central 90% of the counts between 50 and 300 keV have 
been accumulated. 
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MeV. Because this peak is the only one evident at all energies, we chose to run a dedicated 
spectral analysis between T + 9.8 s and T + 10.5 s as described in section H~T1 

We estimated the variability time scale using the full width at half maximum (FWHM) 
of the bright pulse seen around T + 10 s. A combinat ion of exponential functions is used to 
fit the light curve as performed in iNorris et al.l ( 119961 ). The light curve fitting is performed 
for all bright Nal detectors (N6, N7, N8) with a 2 ms time resolution. Two exponential 
functions are used to represent the weak and main bright peaks and include a quadratic 
function to fit the longer timescale variations. As a result, we obtain a FWHM of the main 
peak of 0.15 ± 0.01s for the bright pulse. 



4. Spectral Analysis 



4.1. LAT and GBM spectral fitting 

We performed a time-integrated joint spectral analysis of the LAT and GBM data for 
the prompt phase defined as T + 3.3 s to T + 21.6 s in figure [U For the GBM, we used 
'T ime Tagged Events' (TTE) data from the Nal detectors 6, 7, 8 and BGO detector 1. As 



m 



Abdo et al.l (l2009dl ). background rates and errors are estimated during the prompt phase 



by fitting background regions of the light curve before and after the burst. We derived our 
background estimates using the time intervals [T — 44; T — 8] s and [T + 36; T + 100] s for 
the Nal detectors, and [T — 43; T — 16] s and [T + 43; T + 300] s for the BGO detector. 
For the LAT, we extracted 'transient' class data from an energy -dependent acceptance cone 
around the burst position, as desc ribed in lAbdo et al.l (j2009dl ). and considered front- and 
back-converting events separately ( Atwood et al. 20091 ) . The data files for the analysis were 
prepared using the LAT ScienceTools-v9rl5p2 package, which is available from the Fermi 
Science Support Center (FSSC), and the P6_V3_TRANIENT response functions^ A synthetic 
background was derived for the LAT data using an empirical model of the rates expected for 
the position of the source in the sky and for the position and orientation of the spacecraft 
during the burst interval. 

The joint spect r al fitting of GBM and LAT data was performed using rmfit version 3.2 
(IKaneko et al.ll2006l ; lAbdo et al.ll2009dl ). which estimates the goodness-of-fit in terms of the 
Castor Statistic (C-STAT) to handle correctl y the small number of events at the highest 
energies. The Castor statistic (IDormanl 120031 ) is similar to the Cash statistic (ICashlll979l ) 
except for an offset that is constant for a particular dataset. A global effective area correction 



http : //f ermi . gsf c .nasa. gov/ ssc/ 
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Fig. 1. — GBM and LAT light curves for the gamma-ray emission of GRB 090926A. The 
data from the GBM Nal detectors were divided into soft (8-14.3 keV) and hard (14.3- 
260 keV) bands to reveal similarities between the light curve at the lowest energies and that 
of the LAT data. The fou rth panel shows all LAT events that pass the on-board GAMMA 
filter ( lAtwood et al.ll2009l ). The first four light curves are background-subtracted and are 
shown for 0.1 s time bins. The fifth and sixth panels show LAT data 'transient' class events 
for energies > 100 MeV and > 1 GeV respectively, both using 0.5 s time bins. The vertical 
lines indicate the boundaries of the intervals used for the time-resolved spectral analysis, 
T + (0, 3.3, 9.8, 10.5, 21.6) s. The insets show the counts for each data set, binned using 
these intervals, to illustrate the numbers of counts considered in each spectral fit. 
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Fig. 2. — GBM and LAT light curves for the gamma-ray emission of GRB 090926 A with 
0.05 s time binning for the core of the prompt phase. The vertical dashed lines at T + 9.8 s 
and T + 10.5 s define interval c used in the spectral analysis, see section ED 
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has been applied to the BGO data to match the model normalizations given by the Nal data; 
this correction is consistent with the relative uncertainties in the GBM detector responses. 
Uncorrected, this will normally cause a mismatch between the fitted model rates between 
the two types of detectors where they overlap in energy Once the correction has been 
determined, it is held fixed throughout the calculation, since it reflects an uncertainty in the 
response rather than in the data. In this analysis, the Nal to BGO normalization factor was 
found to be 0.79. For further details on the data e xtra ction and spectral analysis procedures 
see our previous publications lAbdo et al.l ( ]2009dl ) and lAbdo et al.l ( 120101 ) . 



Initially, we fitted a canonical Band function (IBand et al.lll993l ) to the data and then 
found that adding an extra power-law component improved both the fit statistics and residu- 
als. Table [T] summarizes the best-fit parameters and shows that the improvement in C-STAT 
for the (Band+PL) fit over the Band fit alone is 107.3, indicating a firm detection of the 
additional power-law component. The parameters of the Band function are stable, and the 
power-law photon index of the additional component is A = —1.79 ± 0.02. 

In order to better characterize the power-law component at the highest energies, we 
ran a LAT-only data analysis using the unbinned likelihood technique for the full prompt 
phase. The fitted spectrum is shown in figure [3] (black points). The resulting photon index 
is -2.29 ± 0.09, much softer than the -1.79 ± 0.02 index found for the joint GBM/LAT 
analysis. Considering the systematic effects in both analyses, this difference in photon index 
is significant (~ 3a level) and is an indication of the presence of a spectral break. With 
the LAT data alone, we could not find any significant evidence for a deviation from the 
simple power-law shape, probably because of the limited lever arm in energy. Hence, we 
investigated this effect using the joint fits of the GBM and LAT data. 

We fitted the GBM/LAT spectra with the combination of the Band function and a 
power- law model with an exponential cutoff (CUTPL), 



f(E) = B 



E 



cxp 



E 
E~ F 



(1) 



Here B is the normalization in units of photons s 1 cm 2 keV 1 , E P i v is the pivot energy 
fixed at 1 GeV, Ep is the e-folding energy, and A is the power-law photon index. 

The fit results are summarized in table [IJ and the count spectra and residuals are 
shown in figure H] for the best-fit model. The e-folding energy is Ep = 1A1^^ stat. ± 
0.30 syst. GeV, while the power- law photon index below the cutoff energy is A ~ — 1.72+q;q2 stat. ± 
0.01 syst., which is a bit harder than in the (Band+PL) case. The systematic uncertainties 
ha ve been deri v ed usin g the bracketing instrument response functions, as described in detail 
in 



Abdo et al.l (J2009dj). The parameters of the Band function change little from one fit to 
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another. The C-STAT value for this model improves by 40.5 compared to the (Band+PL) 
model, which is significant at the > 4cx level (see the deeper discussion below). We also tried 
to fit the data with a broken power-law model, 



where A; and A^ are the low- and high-energy power-law photon indexes, respectively, E piv 
is the pivot energy fixed at 1 GeV, and -Ebreak is the break energy. However, the significance 
of the fit was close to that found using the (CUTPL) model so that we cannot distinguish 
between the two models. The fit with a broken power-law gave a break energy -Ebreak = 
219t| MeV and a high-energy photon index of X h = -2.471%$. 

One may assess the significance of the spectral cutoff by computing the difference in the 
best-fit C-STAT values for the (Band+PL) and (Band+CUTPL) models. Since C-STAT is 
equal to twice the log-likelihood, this is the standard likelihood ratio test; and conventionally, 
one calculates the significance of a change in log-likelihood using Wilks' theorem. In this 
case, Wilks' theorem states the A(C-STAT) values should be asymptotically distributed as 
X 2 for one degree of freedom. However, certain assumptions are required for the validity of 
this calculation. For the highest reliability, we studied the distribution of A (C-STAT) values 
via simulations, creating 2 x 10 4 random realizations of the null hypothesis (the (Band+PL) 
model with parameters set at the best-fit values) and fit the data for each trial with both 
models. In the resulting distribution of A(C-STAT) values, the largest difference we found 
was 16.7, much smaller than the value of 40.5 for the actual data (see table [1]). We therefore 
place a firm upper-limit on the probability that our fit of the exponential cutoff occurred by 
chance of 5 x 10~ 5 . This corresponds to a Gaussian equivalent significance of 4.05a. 

Our distribution of A (C-STAT) values shows a slight excess over the \ 2 distribution 
at large values indicating that perhaps the asymptotic distribution has not been reached 
for this number of trials. To be conservative, we do not evaluate the significance according 
to the conventional procedure of using the observed A(C-STAT) value of 40.5 and the x 2 
distribution. Unfortunately, the number of simulations that would be required to determine 
the significance of the observed cutoff is prohibitive. Nonetheless, the sizeable gap between 
the largest A(C-STAT) value obtained in the simulations, 16.7, and the observed value of 40.5 
suggests that the significance is much larger than 4cr. For the 4 different sets of instrument 
response functions that we used in our study of the systematic uncertainties, we always found 
A(C-STAT)> 32. The significance of the spectral cutoff will be hereafter quoted as > 4cr. 



Using the fit results for the best model (Band+CUTPL), we estimate a fluence of 
2.07+0.04 x 10" 4 ergcm- 2 (10 keV-10 GeV) from T + 3.3 s to T + 21.6 s. These data 




1 
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give an isotropic en ergy £ 7i j SO = 2.24 ±0.04 x 10 54 erg, comparable to that of GRB 090902B 



f lAbdo et al.ll2009af ) 



We then performed a time-resolved spectral analysis of the prompt phase in the four 
time intervals a, b, c, d. The spectra are shown in figure [5j and the results are summarized in 
table El where the best-fit parameters are given for the statistically preferred model, and the 
C-STAT values are given for the various models. The extra power-law component is found 
to be very significant in intervals c and d, but not at the beginning of the prompt phase in 
intervals a and b. The spectral cutoff is significant at the > 4a level only in the common 
sharp peak (time interval c), where the GeV flux is the highest, but is only marginally 
significant (~ 4cr) in time bin d. 

In time interval b, the improvement in the fit statistics when adding the extra power-law 
component is only A(C-STAT) = 11.6. As a consequence, the parameters of the power-law 
are not very well constrained, yielding a normalization B = 2.9if q 10 -10 photons cm -2 s _1 
keV -1 and a power-law index A = In time interval c, we found the cutoff energy 

to be E F = 0.40+o;oe stat. ± 0.05 syst. GeV (table |2J). Note that we fixed the pivot energy 
at E p[v = 1 MeV for time interval c, since this is the only interval where the extra power- 
law component is dominant over the Band component at very low energies, and setting 
E piv = 1 GeV resulted in very asymmetric and very large uncertainties, especially for the 
normalization B of the extra power-law component. We also tried to fit time interval c with 
a broken power-law model, see equation [2J but again the fit significance was close to that 
of the (Band+CUTPL) model, so that we cannot distinguish between the two models. The 
fit with a broken power law gave a break energy -Ebreak = 264^y5 3 MeV and a photon index 
above -Ebreak of \h = — 3.55I328. In time interval d, the improvement in the fit statistics 
when adding a cutoff to the the extra power-law component is only 17.4 (roughly ~ 4a), 
which is quite high, but not sufficient to claim the presence of an energy cutoff in this 
bin alone. However, as the cutoff is strong in the preceding time interval c, we looked 
at the behavior of the e-folding energy. For time interval d, the e-folding energy is found 
to be E F = 2.21±g;|| GeV, which is much higher than the one found in interval c (the 2-a 
confidence intervals for the cutoff in bins c and d actually exclude each other). This indicates 
a possible time evolution of the high energy cutoff. 



4.2. LAT extended emission 

As the burst was occulted by the Earth from T + 540 s to T + 3000 s, we performed the 
unbinned likelihood analysis using 'transient' class events in the time interval from To + 20 s 
to T + 300 s (a small margin is needed to safely define a circular ROI), and use 'diffuse' 
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class events after To + 3000 s0. 'Transient' class events are treated as in §4.11 In addition, 
for the 'diffuse' class events, we included in the model the standard galactic background 
component, described by the FITS model file gll_iem_v02 . f it, with fixed normalization, 
and the standard isotropic background component, whose spectrum is given in the model 
file isotropic_iem_v02.txt, with the normalization left free. Both model files may be 
downloaded from the FSSC website. 

We divided the LAT data into several time intervals, using intervals a,b,c,d for the 
prompt phase, and modeled the GRB extended emission spectrum as a power-law. For the 
period T + 3000 s - T + 4800 s, the fit resulted in a test statistic of 29.4, corresponding to 
a detection at a ~ 5a level, which is remarkable for a time period ~ 1 hour after the burst. 

Figure [6] shows the flux and photon index versus time. The LAT flux follows a power-law 
with time-dependence (T — T n )~ 1, 69=l=0 03 after T n + 21 .6 s, similar to the behavio r of bursts 
GRB 090510 and GRB 090902B kbdo et all l2009al lcl: kckermann et all boiOal ). Prior to 
To + 21.6 s, the photon index varies significantly with values ranging from —2.5 to —1.7. By 
contrast, after To + 21.6 s, the photon index is almost constant with values in the range —1.5 
to —1.9. The soft spectral index in time interval c is consistent with the spectral break of 
the extra component described in section 14. 1[ and the gradual hardening from time bin d is 
consistent with its disappearance. 



1 See lAtwood et al.l ([20091 ) for the definitions and recommended usage of the LAT event classes. 
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Table 1: Summary of GBM/LAT joint spectral fitting between T + 3.3 s and T + 21.6 s. 
The flux range covered by both instruments is 10 keV-10 GeV. 



Fitting model 


Band 


Band+PL 


Band+CUTPL 


Band function 








A (7 cm -2 s _1 keV -1 ) 


0.176 ± 0.002 


0.173 ± 0.003 


1 70 +0 - 001 

u - inj -0.004 


£ P cak (keV) 


249 ± 3 


256 ± 4 


259 ±\ 


a (index 1) 


-0.71± 0.01 


-0.62± 0.03 


-0 64+ - 02 

u -°+-0.09 


/3 (index 2) 


-2.30± 0.01 


_ 2 cq+0.04 


O fi O + 0.02 

z - uo -0.12 


Power-law 








B (1(T 10 7 cm" 2 s" 1 keV- 1 ) 




1 7 +0.35 
01 ' -0.33 


5 80 + - 81 
°- ou -0.60 


A (index) 




-1.79± 0.02 


1 7 9 +0.10 
L - ' ^-0.02 






1 GeV (fixed) 


1 GeV (fixed) 


High- energy cutoff 








E F (GeV) 






1 41 +0-22 
-0.42 


Flux (7 cm~ 2 s _1 ) 


42.2±0.1 


43.5±0.3 


43.3±0.2 


Flux (10~ 5 erg cm" 2 s _1 ) 


1.18±0.01 


1.15±0.02 


1.13±0.02 


C-STAT / DOF 
A(C-STAT)f 


1395.1 / 579 


1287.8 / 577 
107.3 


1247.3 / 576 
40.5 



f with respect to the preceding model (column). 
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Table 2: Summary of GBM/LAT joint spectral fitting by best model in 4 time intervals. The 
flux range covered by both instruments is 10 keV-10 GeV. 



Time interval from To (s) 


(a) 0.0-3.3 


(b) 3.3-9.8 


(c) 9.8-10.5 


(d) 10.5-21.6 


Prefered Model 


Band 


Band 


Band+CUTPL 


Band+PL 


Band function 










A (7 cm" 2 s- 1 kcV- 1 ) 0.146 ± 0.004 


0.302 ± 0.004 


n Qoc +0.064 
U.OOO _ .012 


0.100 ± 0.005 


T-peak 


338 ± 10 


288 ± 4 


209 ±5 e 


186 ± 6 


a (index 1) 


-0.42 ±0.03 


-0.55 ±0.01 


_n 5Q+ - 39 

u - oy -0.06 


-0 70+ - 07 

u - ' u -0.06 


f3 (index 2) 


z -°^-0.09 


-2.46± 0.02 


-3 69+ 1 81 


-2 80+ 13 


Power-law 










B (10- 10 7 cm" 2 s- 1 kcV- 1 ) 






7 56 t +2.25 
'• 0D -0.50 


A7 +0.38 
°- u ' -0.36 


A (index) 






1 71+O.O2 
L - 1 1 -0.05 


-1.79 ±0.03 


R ■ 






1 MeV (fixed) 


1 GeV (fixed) 


High- energy cutoff 










E F (GeV) 






40 

u -^ u -0.06 




Flux (7 cm -2 s _1 ) 


31.4 ±0.2 


66.4 ± 0.3 


109.0 ±1.4 


25.5 ±0.2 


Flux (lO -6 erg cm -2 s _1 ) 


9.96 ±0.41 


18.9 ±0.30 


29.22 ±1.60 


5.83 ±0.30 


C-STAT / DOF 










Band 


622.4 / 579 


944.2 / 579 


655.9 / 579 


1033.8 / 579 


Band+PL 


624.3 / 577 


932.6 / 577 


598.7 / 577 


950.6 / 577 


Band+CUTPL 


618.8 / 576 


928.3 / 576 


574.2 / 576 


933.2 / 576 


A (C-STAT) 










Band^ (Band+PL) 


-1.9 


11.6 


57.2 


83.2 


(Band+PL)^(Band+CUTPL) 


5.5 


4.3 


24.5 


17.4 


^ As E piv = 1 MeV ; B has the unit 


of 1Q- 4 7 cm' 


2 s- 1 keV' 1 
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Fig. 3. — vF v spectrum of the data points from the LAT-only unbinned likelihood analysis 
of GRB 090926A between T + 3.3 s and T + 21.6 s. Black dashed and solid lines show the 
best-fit power-law model and ±1 a error contours, derived from the covariance matrix of the 
fit. 
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Fig. 4. — Joint spectral fitting of GBM and LAT data between To + 3.3 s and To + 21.6 s. 
The top panel shows the count spectra and best-fit (Band+CUTPL) model (histograms). 
The lower panel shows the residual of the spectral fitting. 
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Fig. 5. — Top: The best-fit (Band+CUTPL) model for the time-integrated data plotted as 
a vF v spectrum. The two components are plotted separately as the dashed lines, and the 
sum is plotted as the heavy line. The ±1 a error contours derived from the errors on the 
fit parameters are also shown. Bottom. The vF v model spectra (and ±1 a error contours) 
plotted for each of the time bins considered in the time-resolved spectroscopy. 




Fig. 6. — Time variation of the LAT flux (red cross) and photon index (blue filled circle) for 
the extended emission of GRB 090926A . After the end of the prompt emission at T + 21.6 s, 
the flux decays following a power-law with index —1.69 ± 0.03 (red solid line). 
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5. Discussion and Interpretation 



5.1. Prompt Emission Phase 

The Fermi observations show that GRB 090926A clearly has an extra high-energy com- 
ponent in addition to the Band component in the time-integrated as well as in the time- 
resolved spectra. This is the third case of a LAT detection of an extra spectral com ponent, 
after GRB 090510 and GRB 090902B (kckermann et alihoiOal : kbdo et aDl2009ah . That 
is, we have such detec tions in 3 out of the 4 brightest LAT GRBs (except GRB 080916C, 



sec 



Abdo et al.ll2009bl ). Since we require a confidence level of > 5<r to claim a detection, we 



can unambiguously identify it only in the brightest LAT GRBs, which suggests that such a 
component may be intrinsically very common in GRBs. 

The behavior in time interval c is remarkable, as the light curve shows a clear spike- 
like structure in the LAT energy range that is dominated by the extra spectral component. 
The spike in all the energy ranges has the peak times synchronized within 50 ms and shows 
similar pulse widths (see Figure [2]), which indicates that the origins of the lowest and highest 
energy emission components as well as the Band component are related, i.e., they could either 
arise from the same physical region and possibly also the same spectral component, or are 
otherwise directly physically linked, such as photons generated in one emission region being 
scattered by electrons in the other emission region (see specific models discussed below). 

The delayed o nset of the LAT emission is common to almost all the LAT GRBs (except 
GRB 090217A, see lAckermann et al.ll2010bl ). The delay may arise from the following four 
effects: (1) a flux increase of the Band component, (2) a hardening of the Band component 
(i.e., increase of the peak energy -E pea k and/or high-energy spectral index 0), (3) a flux 
increase of the extra component, or (4) an increase of the cutoff energy in the spectrum. 
However, effect (4) does not seem to be a major effect for LAT GRBs so far since there is no 
clear sign of a high-energy cutoff o r steepening in the spectra before the LAT onsets (see also 
the discussion on GRB 080825C in kbdo et al1l2009dh . The LAT detection of GRB 090926 A 
starts from time interval b, and the clear emergence of the extra component occurs even later, 
10 s after the onset of the Band component, so that this delay is likely to be due to the 
combination of (1) and (2). However, while no significant extra component is detected in 
time intervals a and b, it may still be present with a lower cutoff energy that falls under 
or close to the Band function model component, so that we may not exclude contributions 
from effects (3) or (4) to the delayed onset. 

There are several theoretical models for the origin of the extra spectral component. 
The delayed extra co mponent could be emitt e d from a forwar d shock that propagates into 
the external medium (IMeszaros fc Reed 1 199 71 ; ISari et al.lll998l ). while the Band component 
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is thought to have a separate origin. The delay timescale of the extra spectral component 
would correspond to the time needed for the forward shock to sweep up material and brighten 
( IKumar &: Barniol Duranll2009l ; iGhisellini et al.ll2010t lRazzaqudl2010f ) . The rapid variability 
observed in GRB 090926A is contrary to expectations from an external shock model, unless it 
is produced by emission from a small portion of the blast wave within the Doppler beaming 
cone. This could occur, for instance, if the external medium is clumpy on length scale 
« T / cAT/(l + z) ~ 10 12 (r//10 3 )(AT/0.2 s) cm, where T f is the Lorentz factor of the 
forward shock and AT is the pulse duration ( iDermer fc Mitmanlll999l ; lDermerll2008l ). This is 
based on interactions between a very thin shell, prior to the onset of the self-similar expansion 
phase, and an external medium with very small scale clumps. If the extra component is 
synchrotron emission from the forward shock, then the synchronization of the pulse peak 
times of the Band and extra component requires an explanation. One possibility is that the 
extra component arises from inverse Compton (IC) scattering of the radiation of the Band 
component by the high-energy electrons in the forward shock. 

As for internal emission models, in which both spectral components arise within the 
ejecta, the extra component can be produced by IC scattering by energetic leptons or via 
hadronic processes. In either case, the time of the peak of the extra component would lag 
relative to the Band component in the same emission episode, although the time lag can be 
limited by a timescale comparable to the pulse duration, which would still be consistent with 
the observed synchronization of the two components. A simple leptonic model could comprise 
synchrotron plus synchrotron self-Compton ( SSC) emission. Under t he physical conditions 
typically assumed in the internal shock model (IRees &: Meszaroslll994j ). all electrons emitting 
synchrotron emission cool on a timescale much shorter than the dynamical time (i.e., the 
electrons are in the fast cooling regime), so that the photon index a below -Ep ea k should 
be —1.5, which is not consistent with the results from our fits, —0.7 < a < —0.4. The 
synchrotron plus SSC model would need to overcome this problem . In the photospheric 
emission model of the Band component (e.g., IMeszaros fc Reed 120001 ). the extra component 
could b e the IC of the photospheric emission by electrons in the dissipation region at large 
radius (lBeloborodovll2009t iGao et al.ll2009t iRyde et al.ll201Q[ ). In this model, the low energy 
excess seen in interval c could be synchrotron emission of the electrons and the electron- 
positron pairs created by the cascade process, and the delay timesca le of the extra com ponent 
could be explained by the evolution of the jet physical conditions (IToma et al.ll2010l ). 



Hadr onic processes, such as a photopion- induced pair cascade or p roton/ion synchrotron 
emission ( lAsano et al.l 120091 ; iRazzaque et al.l I201CH ; IWang et al.l 120091 ) can make a spectral 
component that is distinct from that which is commonly observed during the prompt phase. 
Large values of the bulk Lorentz factor of the emitting region, T > 10 3 , imply large energy 
requirements for significant high-energy emission in either photo-hadronic or proton/ion 
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synchrotron models. The bulk Lorentz factor that is inferred by the argument of the pair 
absorption opacity for this burst, however, is relatively low (r ~ 200-700, see below). Thus 
for proton/ion synchrotron models, whi ch require a heavi l y magnetically loaded shocked jet, 
the total energy requirements oc T 16 / 3 ( IWang et al.ll2009t iRazzaque et al.ll2010l ) are smaller 
by a factor of ~ 7-5000 compared to the case of T ~ 10 3 , much improving the viability of 
such models. The lower allowed values of T also reduce the energy requireme nts in photo 



hadronic models, where non-thermal protons usua lly dominate the bulk energy (lAsano et al. 



20091 ; IWang et al.l 120091 ; lAckermann et al.l l2010a[ ). Under the assumption that the photon 
field is homogeneous and steady in the emitting region, the efficiency of photo-hadronic inter- 
actions by high-energy protons at the peak of the Band spectrum component in this burst is 
estimated to be a few percent for T ~ 10 3 . This efficiency scales as T -4 , which implies signif- 
icant reduction of the total energy requirements for the lower allowed values of T. However, 
r < 700 is indicated under the assumption that the photon field is inhomogeneous and time- 
dependent, as will be discussed below, for which it is not so clear how the photo-hadronic 
i ntera c tion efficiency ( as we ll as the neutrino production efficiency, see iMurase fc Nagataki 
tood ): bermer et aD t00l\ ): iRazzaquel hood \) depends on T. 



Another remarkable aspect of this burst is the spectral break (or cutoff) of the extra 
component that has been measured in the time-integrated spectrum of the prompt emission 
and for time interval c with a high significance (> 4 a; see sections |4.1[ and tables [U E]). This 
cutoff may be due to pair production (77 — > e + e~) within the emitting region, although we 
cannot rule out the possibility that there is an intrinsic spectral break related to the energy 
distribution of the emitting particles or the emission mechanism (e.g., IC scattering in the 
Klein-Nishina regime). Absorption by the extragalactic background light (EBL) cannot 
cause this spectral feature since the opacity at the observe d break energy for the redshift of 
GRB 090926A is very small for practically all EBL models (IFinke et al.ll2010l . and references 
therein). We focus on the spectral feature in time interval c to constrain the physical 
properties of the emitting region by introducing t he critical photon energy E n at which the 
pair production opaci ty is unity, r 77 (i? c ) 
Lithwick fc SarilboOlh . 



1 (e.g.. Krolik fc Pierlll99ll ; iFenimore et al.lll993 



In order to characterize the spectral break, we have fit the data with a model that 
consists of an extra power-law component modified by absorption due to pair production, 
but this spectral model is not unique. Although the instantaneous emission from a thin 
shell exhibits a photon spectrum like / oc E x exp(— r 77 (£■)), the shape of the time-int egrated 
spectrum of a single pul se may depend on the details of the emission mechanism (jBaring 
20061 ; iGranot et al.ll2008l ). For example, the simple model of an emitting slab leads to / oc 
E x [l — exp(— T yy (E))]/r^(E), which is a smoothly broken power-law spectrum since r^(E) 
is a power-law function of E when the intrinsic emission spectrum is a power-law function 
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(see below). A fully time-dependent and self-consistent semi-analytic calculation featuring 
emission from a very thin spherical shell over a finite range of radii (IGranot et al.l 120081 ) 
would also lead to a smooth break to a steeper power-law in the time-integrated spectrum of 
a single pulse. In section 4.2, we have fitted the observed extra spectral component for time 
interval c by two empirical functions: a power-law with exponential cutoff (Eq. 1) and a 
broken power-law model (Eq. 2). However, the photon counts are not enou g h to d istinguish 
between the two models. In the time-dependent model of IGranot et al.l (120081 ). photons 
above the spectral break energy are expected to arrive predominantly near the onset of the 
spike. This signature may afford a more unambiguous indicator of the intrinsic opacity to 
pair production. Such an opportunity is also unavailable for GRB 090926A, however, due to 
insufficient photon counts above the spectral break. Here we only use the result of the former 
model (Eq. 1), and consider the e-folding energy Ep from the fit to be good approximation 
of E r . 



In order to derive the pair absorption function t^(E), we first consider a simple model 
in which the photon field in the emitting region is unif orm, isotropic, and time-independent 
in the comoving frame (see the supporting material for lAbdo et al.ll2009bl ). We assume that 
the opacity at the photon energy around E c is dominated by the extra power-law component 
itself, instead of the Band component. This assumption is justified for the observed spectrum 
in time interval c, as shown below. Let us define the observed photon number spectrum of 
the extra component for one pulse, below the break energy, as f(E) = f (E piv )(E / E piv ) x in 
units of photons cm" " 2 keV -1 . The energy distribution of the photons in the comoving frame 
of the emitting region is written as 



<{E>) 



R 



piv J 



'l + z) 3 W \E' 



E' 



piv 



(3) 



where the quantities with a prime are measured in the comoving frame, ~ 5.17 x 10 cm 
is the luminosity distance of the source, and R, T, and W are the distance from the 
central engine, the bulk Lorentz factor, and the comoving radial width of the emitting 
region, respectively. Photons with energy E' = E' c > m e c 2 annihilate mainly with tar- 
get photons with energy E' ann ~ 2m 2 c 4 / ' E' c . Then the optical depth is of the order of 
T yy (E' c ) ~ Q.laTE'^n'JE'g^W' , where the pair production cross section is approximated to 
be 0.1 times the Thomson cross section op- More accurately, we have 
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m 2 c 4 



-A-l 



F(X) 



(4) 



f lAbdo et al.l l2009bh where F(X) « 0.597(-A)~ 2 ' 30 for -2.9 < A < -1.0. The relation 
R ~ T 2 cAT/(l + z) is valid for a large class of emission mechanisms, where AT is the 
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variability time. Then we obtain 

i 

/ TP TP \ ~ X ~ l \ 2(1 ~ A) 

(5) 



i_ 



We can estimate Y and i? from the values of E piv , f(E piv ), A, and E c , which are inferred from 
the observed spectrum. We have estimated the variability timescale of the Band component 
to be AT ~ 0.15 ± 0.01 s from the analysis of the GBM emission in section 3. From 
the synchronization of the peak times and the similar pulse durations in all energy ranges, 
we may assume that the variability timescale of the extra component is the same. The 
power-law with exponential cutoff model with E piY = 1 MeV results in B = 7. 56^50 x 
10" 4 photons cm" 2 keV~\ A = -1.71±°;° 2 , and E c = E F = 400+^° MeV. We may estimate 
the normalization of the spectral fluence over the variability time AT around the spike by 
f(E piv ) ~ 2BAT, where B is the normalization of the time-averaged spectral flux over 
interval c. Then we obtain V ~ 720 ± 76, where the error on V has been calculated by the 
statistical errors of the parameters X,Ep, and AT as well as the error of f(E piv ) taken as 
±5 AT/2. This error on T is much smaller than the uncertainty of T due to modeling the 
space inhomogeneity and the time dependence of the target photon field discussed below. The 
energy of the main target photons for the photons with E c is -E aim ~ 2r 2 m 2 c 4 /[(l + z) 2 E c ] ~ 
70 (r/700) 2 (Tv c /400 MeV) -1 MeV, while the extra power-law component is dominant above 
~ 1 MeV. Thus, our assumption that the target photons for the photons with energy E c are 
from the extra component is justified. In this model, the emission radius is estimated to be 
R ~ r 2 cAT/(l + z)~7x 10 14 (r/700) 2 (AT/0.15 s) cm. 



A fully time-dependent and self-consistent semi-analytic model by iGranot et ajj (120081 ) 
results in a significant reduction in r 77 and in the inferred value of T by a factor of ~ 3 
compared to simpler models like the above calculation, i.e., T ~ 220 for this burst. Under 
the typical physical conditions for the shock emitting the bright 7-rays, electrons are in 
the fast cooling regime, so that most of the radiation is em itted within a very thin layer 
behind the expanding shock front (e.g., IGranot et al.ll2000l ). The reduction in r 77 occurs 



mainly since the high-energy photons are emitted from a very thin cooling layer, so that 
those that are emitted from angles < 1/T relative to the line of sight immediately propagate 
ahead of the shock front, and can therefore potentially pair produce only with photons that 
propagate at fairly small angles, #12, relative to their own direction. The small values of 
6*12 suppress the interaction rate (r 77 oc 1 — cos6 , 12 ), and increase the threshold energy for 
pair production, EiE 2 (l + z) 2 > 2(m e c 2 ) 2 /(l — cos#i 2 ), where E\ and E 2 are measured at 
Earth. The time dependence also reduces the time-averaged opacity over a single spike in the 
light curve, since the opacity is initially very low and gradually increases as the photon field 
builds-up, approaching its quasi-steady state value on the dynamical time (by which time 
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the emission episode leading to the spike in the light curve is typically over). Furthermore, 
in this model both the photon density and the typical value of 9 12 decrease with radius along 
the trajectory of a test photon, further reducing r 77 . For our discussion below we adopt the 
value of T ~ 200-700, intermediate between the values inferred by the previous simple model 
and the time- dependent model for a very thin cooling layer. The motivation for this is that 
in some of the models discussed above the high-energy photons are expected to be emitted 
from the bulk of the shocked region, rather than from a thin cooling layer behind the shock 
front, in which case such an intermediate value of the opacity might be expected. 

The spectrum for time interval d also has an extra power-law component. It is much 
dimmer than that for time interval c, but its spectral index is similar. The break feature is 
marginally significant. A straightforward interpretation of this behavior is that the critical 
energy E c is larger than that for interval c. However, it is also possible that a different 
emission component, which is responsible for the LAT temporally extended emission at 
T > 20 s, contributes to the high-energy emission in interval d. This may hide a possible 
spectral break of the extra component. 

We emphasize that this burst is the first GRB that exhibits a spectral break that can 
be used to estimate (as opposed to bound) its bulk Lorentz factor, presuming that this 
feature is due to pair production attenuation. Other LAT GRBs do not show any clear high- 
energy spectral brea ks (a circumst ance evinced in EGRET bursts with lower count statistics; 
see the overview in iBaringi 120061 ). The lower limits of the bulk Lorentz factors for those 
GRBs have been derived by r 77 (E h ) < 1, where E h is the highest photon energy detected; 
r > 900, T > 1200, and Y > 1000, for GRB 080916C, GRB 090510 and GRB 090902B, 
re spectively, using the simple model described above. In the time-dependent thin-shell model 
of iGranot et al.l (120081 ) all of these lower limits would be lower by about a factor of ~ 3. 
Thus, the inferred T ~ 200-700 of GRB 090926A is smaller than the lower limits for other 
LAT GRBs. On the other hand, it is consistent with the constraints on T for other GRBs 
put by different methods: the observed broad flux peaks of some optical afterglows, typically 
~ 10 2 -10 3 s after the burst triggers, appear to signify the afterglow onset, and thus indicate 



r-iu s alter tne burst triggers, ap p 
a few x 100 iMolinari et alJl2007t K 



lOates et al.ll2009l ); the possible thermal emission in 
radii of the iets, which indicate T ~ 300-700 



(Pe'er et al. 


2007; 


Rvde et al. 


2010) 



widely distributed over a range > 10 2 . 
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5.2. Temporally Extended Emission 



The emission in the LAT energy range lasts much longer than that in the GBM energy 
range, and the flux in the LAT energy range shows a power-law decay, F v oc (T — To) a 
with a = —1.69 ± 0.03. This behavior is similar to that seen in other LAT GRBs, which 
have decay indice s a = -1.2 ± 0.2 for GRB 080916C JAbdo et alil2009bh . a = -1.38 ± 0.07 
for GR B 090510 (foe Pasquale et alJboioh . and a « -1.5 for GRB 090902B kbdo et al. 



2009 aj). These may be explained as synchrotron emi s sion from the external f orward shock 
rtKumar k Barniol Duranlbood : iGhisellini etaJboiot be Pasquale et alJboioh . 



For GRB 090926A, this interpretation seems consistent with the nearly constant spectral 
index at > 20 s (Figure 6). In this scenario, the starting time of the self-similar phase of the 
forward shock should be < Td ur — 13 s. This means th at the ejecta is in the thick shell regim e 
or in the borderline of the thick and thin shell regimes (ISarilll997t iKobayashi fc Zhanp2003l ). 
which corresponds to T dec = [{^-k)£ ktim /{2^A^(^ k T^ k ^^ k \l + ^ < T dur ~ 13 s, 
where iso is the isotropic-equivalent kinetic energy of the ejecta producing the external 



shock, and the external density profile is defined as n, 



A Pxi R k . This relation allows 



us to put a lower limit on the initial Lorentz factor of the ejecta. For the uniform density 



case (A; = 0), V > 750 (£^ iso /2 x 10 54 erg) 1 /8( r/7 /0.2)- 1 / 8 (n ext /l cm- 



medium case (k = 2), T > 290 {£, 



7, is. 



D /2 x 10 54 erg) 1 / 4 (r/ 7 /0.2)- 1 / 4 (A cxt /5 x 10 11 g cm" 1 



and for the wind 

1/4 



where we define ^ 7 = £y t iso/£k,iso> an d A ext ~ 5 x 10 11 g cm 1 corresponds to the value 

1O~ 5 M yr" 1 and 



of typical Wolf-Rayet stars in o ur Galaxy, which h ave mass loss rates 
wind velocities ~ 1000 km s" 1 . ICenko et al.l ( 120101 ) analyzed the observed late optical and 



X-ray afterglows of GRB 090926A and showed that they can be explained by synchrotron 
emission from the adiabatic forward shock propagating into the wind medium with A cxt ~ 
3.4 x 10 11 g cm -1 . This could provide T > 320 (r^/0.2) -1 / 4 , which may be consistent with 
T ~ 200-700 of this burst inferred by the pair attenuation opacity argument. Note that the 
estimate T ~ 200-700 is only for the bulk Lorentz factor of the shell emitting the spike in 
time interval c, while the lower limit T > 320(?7 7 /0.2)~ 1 / 4 is relevant for the mean or typical 
value, weighed over the energy in the whole outflow. 



6. Conclusions 

GRB 090926A is one of the brightest long bursts detected by the GBM and LAT instru- 
ments on Fermi with high energy events up to ~ 20 GeV. As in other bursts (GRB 090510, 
GRB 090902B), this burst shows an extra hard component in its integrated spectrum, but 
for the first time we significantly detect a spectral break around 1.4 GeV. The time-resolved 
spectral analysis shows that the extra component significantly dominates the emission in 
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the high (> 1 MeV) energy range at the time of the narrow pulse which is simultaneously 
observed by LAT and GBM. At earlier times, the spectrum is described by a standard Band 
model while at later times the extra component is significant, but a spectral break feature is 
only marginally significant. Correlation between the lowest and highest energy light curves 
implies that the origins of the Band component and the extra power-law component are 
related around the time of the sharp pulse. 

The ~ 3.3 s delay of the LAT emission onset can be explained as the overall flux 
increase and the spectral hardening of the Band component, since the clear emergence of the 
extra component occurs only at a later time. However we may not exclude a contribution 
from the extra component in the early times, whose flux is intrinsically just below the 
Band component or suppressed by a lower spectral cut-off. The high temporal variability 
of the extra component and the correlation of the Band and extra components put strong 
constraints on the external shock scenario: the external medium needs to be highly clumpy, 
and the emission mechanisms of the two components should be related. 

From the spectral break we have computed the bulk Lorentz factor of the emitting 
shell and find a range of T ~ 200-700, depending on the assumption of the homogeneity 
and time-dependence of the photon field, as well as on the assumption that the cutoff is 
due to the pair production attenuation. Even if we cannot distinguish between leptonic 
and hadronic emission for the extra component, we note that such a moderate Lorentz 
factor could alleviate the problem of the energy budget in hadronic emission models, as 
for GRB 090510. Comparison of this estimate with the large lower limits for other LAT 
GRBs and the estimates for GRBs that occurred before Fermi may imply that the bulk 
Lorentz factors of GRBs are widely distributed over a range of values > 100. In addition, 
the early deceleration of the forward shock inferred by the LAT temporal extended emission 
and the dens i ty of the external medium inferred by the late optical and X-ray afterglows 



( iCenko et al.l 120101 ) can put a lower limit on the bulk Lorentz factor of the entire shell 
just before the deceleration, T > 290(?7 7 /0.2)~ 1//4 , which is consistent with the estimate of 
T ~ 200-700 for the region that corresponds to the emission around the spike. 

Further LAT detections of bright GRBs will enable us to observe other bright extra 
components and constrain their origins and spectral breaks and their relation to excesses 
below 20 keV and the temporally extended emission. 
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